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ABSTRACT 

Context. BP Cru is a well known high-mass X-ray binary composed of a late B hypergiant (Wray 977) and a neutron star, also 
observed as the X-ray pulsar GX301-2. No information about emission from BP Cru in other bands than X-rays and optical has been 
reported to date in the literature, though massive X-ray binaries containing black holes can have radio emission from a jet. 
Aims. In order to assess the presence of a radio jet, we searched for radio emission towards BP Cru using the Australia Compact Array 
Telescope during a survey for radio emission from Be/X-ray transients. 

Methods. We probed the 41.5 d orbit of BP Cru with the Australia Telescope Compact Array not only close to periastron but also 
close to apastron. 

Results. BP Cru was clearly detected in our data on 4, possibly 6, of 12 occasions at 4.8 and 8.6 GHz. Our data suggest that the 
spectral index of the radio emission is modulated either by the X-ray flux or the orbital phase of the system. 

Conclusions. We propose that the radio emission of BP Cru probably arises from two components: a persistent component, coming 
from the mass donor Wray 977, and a periodic component connected to the accretion onto the neutron star, possibly coming from a 
(weak and short lived) jet. 
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1. Introduction 

The majority of High-Mass X-ray Binaries (HMXB) are classi- 
fied as Be/X-ray transients, i.e. the compact object in these sys- 
tems is the companion of a Be-star Some Be/X-ray transients 
display X-ray outbursts at most of periastron passages as the re- 
sult of mass transfer from the Be-star. Another important group 
of HMXBs contains systems were the mass donor is a OB super- 
or hypergiant, that is losing mass through a massive stellar wind, 
from which the compact object is accreting. 

The wind accreting X-ray pulsar GX301-2 (also known as 
4U 1223-62) is part of the well known high-mass X-ray binary 
BP Cru. This system is peculiar in several aspects (Table [1]): it 
contains the most massive star known in these systems and one 
among the slowest pulsars. 

The orbital period of BP Cru is Porh = 41.5 d, and the X- 
ray fl ares occur about two days before periastron dWatson et alj 
|1982|) . In addition BP Cru shows a weaker X-ray flare close to 
apastron. The fact that the strong X-ray outbursts occur close to 
periastron makes it similar to the Be/X-ray transients. Models 
of fast and dense stellar winds from the primary have difficulties 
in reproducing the coiTelations of the weak and strong X-ray 
outbursts with the orbital phases of BP Cru. Among the alterna- 
tives are the existence of a circumstellar disc around Wray 977 
through which the neutron star passes at every peri astron, an en- 
hance ment of the stellar wind close to periastron dPravdo et alj 
Il995h . as wen as a spiral-type tidal stream of gas trailing the 
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mass donor's rotation jLeahv & Kostkall2008l) . that can produce 
a transient accretion disc around the neutron star, whic h explains 
the e pisodes of rapid spin-up of the X-ray pulsar (iKoh et al.l 
Il997h . 

While Be/X-ray transients are easily observed in the X-rays 
and some of them have known optical counterparts, very few 
radio observations of these objects have been re ported in the 
literatu re. The radio surveys o f X-ray binaries bv iDuIdig et al.l 
d 19791) and iNelson & Spenceii (Il988l) report only upper limits 
for 7 Be/X-ray transients, and one positive detection o f 14mJy 
at 2 cm (April 1978) for A 1 118-61 (iDuIdig et aI.|[T979l) . For this 
source they found only upper limits on the 3 mJy level at 7 other 
epochs. The orbital phases of the system during these observa- 
tions are unclear. The main problem for the observation of these 
systems is the transient nature of the outbursts, and in some cases 
the coarse X-ray positions. 

Some other X-ray binaries (as e.g. microquasars), have rel- 
ativistic jets that produce synchrotron emission at radio wave- 
lengths. In general microquasars are black hole candidates, but 
in a few cases the compact object might rather be a neutron star. 
The neutron star systems that are emitting radio waves tend to 
have high luminosities and weakly magnetised neutron stars. On 
the other hand. X-ray pulsars that have strongly magnetised neu- 
tron stars but are less luminous, are as a rule not detected at radio 
wavelengths. 

Outbursts in microquasars are characterised by an X-ray out- 
burst close to periastron followed by the formation of a radio jet 
at a later time, when the X-ray spectrum indicates that the inner 
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Table 1 : Characteristics of the BP Cru system. 



Table 2: Details of the observations presented in this paper. 



BP Cru 



RA(J2000) 12:26:37.6 
Orbital period 



DEC(J2000) -62:46:14 
41.5d 



Wray 977" 



Mass 


39 -68 Mo 


Radius 


62 Ro 




18100K 


Luminosity 


5 X lO^Lo 


Mass-loss rate 


10-^ Moyr-i 


Wind velocity 


305 km s"' 


GX301-2 


Mass" 


1.85 Mo 



Spin period'' 

Surface magnetic field* 



696 s 
4x lO'^T 



" iKaper et alj|200g '' iKrevkenbohm et alj|2004l ^ IWhiteet all 19761 



part of the accretion disc has vanished (iHarmon et al.ll99"5l) . The 
radio outburst can therefore be delayed by up to some 10 days 
relative to the X-ray outburst. 

With this in mind, we devised an observational programme 
to search for radio emission from Be/X-ray transients and other 
HMXBs that have semi-regular X-ray outbursts close to perias- 
tron. Our strategy was to distribute "snap-shot" observations to 
cover the orbital phases between periastron and apastron. 

In this letter we present the first results of our observational 
campaign, the detection of radio continuum emission in the cen- 
timetre band at two frequencies (4.8 and 8.6 GHz) from the 
BP Cru system. 

2. Observations and data reduction 

Observations of BP Cru were done at 12 epochs for a total of 
approximately 20 hours of observations, 13 of which on target. 
The observation epochs were spread between November 2008 
and February 2009, covering almost four complete orbits. Every 
epoch consisted of about 1.5 hours of data recording, and this 
time was divided between flux and phase calibrators and the tar- 
get. The original strategy was to spread our observations over 
different orbital phases, from shortly before periastron to close 
to apastron. Our proposed schedule was then folded with the 
telescope availability to give the final schedule summarised in 
Tableland graphically presented in Fig.[T] 

We used the Australia Telescope Compact Array (ATCA) 
in the 4.8 and 8.6 GHz continuum bands in two linear polari- 
sations. Data were recorded at the largest bandwidth, 128 MHz, 
spread over 33 channels per polarisation and frequency. All four 
polarisation products were recorded (XX, YY, XY, YX) in or- 
der to retrieve information about the level of polarisation of a 
detection. The time in every observation block was divided be- 
tween a phase calibrator (1148-671), and the target, alternating 
2 minutes scans on the calibrator with 6-8 minutes on the target. 
Additionally, at the beginning of each run, a bandpass calibra- 
tor was observed for some 5-10 minutes. Since the project was 
a detection experiment, no requirements were made as to which 
configuration ATCA had to be in at the different epochs. This had 
the advantage of greatly improving the scheduling flexibility but 
the disadvantage of having to reduce data in all kinds of config- 
urations. In this respect, data at epochs 4, 5, 6, 11 and 12 (see 



Epoch 


Date (UT)" 


T-on [min]* 


Phase^ 


Array"* 


1 


2008-11-01, 2212 


36 


-7 


6A 


2 


2008-11-05, 1708 


54 


-3 


6A 


3 


2008-11-09, 0443 


63 


+ 1 


6A 


4 


2008-11-18, 0314 


45 


+ 10 


EW367 


5 


2008-11-24, 2055 


45 


+ 16 


EW367 


6 


2008-12-16, 1936 


54 


-4 


750B 


7 


2008-12-21, 1840 


54 


+ 1 


6C 


8 


2009-01-05, 1719 


153 


+ 16 


6C 


9 


2009-01-26, 1603 


63 


-4 


1.5C 


10 


2009-01-30, 1415 


72 


+0 


1.5C 


11 


2009-02-08, 1936 


72 


+9 


EW352 


12 


2009-02-15, 1122 


54 


+ 16 


EW352 



" Times indicate the beginning of every observation 

* Effective time spent on-source 
Days before (-) or after (+) periastron. Peria stron passage occ urred 
at 2008-11-08, 2008-12-20, and 2009-01-30 (see lKoh et alll997h . 

'' ATCA configuration for the observations. When including all an- 
tennae, the Mv-coverage from the '6' -configurations is regular, from 
the ' 1.5'-configurations slightly bimodal and from the EW- and 750- 
configurations (compact configurations) it is strongly bimodal. 



Nov 2008 



V .. 



Jan 2009 Feb 2009 

Date 



Fig. 1: X-ray lightcurve (Swift/BAT, 15 - 150 keV) of BP Cru 
in the period of our observing campaign (red dots). Vertical 
black lines indicate the epochs of our radio observations, and 
the dashed blue lines the periastron passages. Green triangles 
peak up indicate a positive spectral index, magenta triangles 
peak down indicate negative spectral index. 



Table |2]) turned out to be less useful for our program because 
they were recorded in the most compact ATCA configurations 
(see Table|2|i and are hence heavily biased by the Galactic emis- 
sion, as BP Cru lies in the plane of the Galaxy. The 1 cr rms level 
in those maps was on average ^ 0.3 mJy (i.e. a detection limit of 
~0.9mJy) at both frequencies. To prevent loss of sensitivity on 
target due to possible artefacts, the ATCA array was pointed ±10 
arcsec away from BP Cru. Finally, most of the observations were 
done in remote mode, making use of the exceptionally functional 
software available for ATCA. 

Bandpass, flux and phase calibration were performed using 
standard procedures in AIPS ' with a ParselTongue^ interface to 
allow quick and equal reduction of many observational epochs, 
as well as in MIRIAD^. The bandpass calibrator was 0823-500 



http ://w w w. aips.nrao.edu/ 

http://www.radionet-eu.org/rnwiki/ParselTongue 

http://www.atnf.csiro.au/computing/software/miriad/ 
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Table 3: Details of the detections and upper limits of the non- 
detections of BP Cru ordered in orbital phase. 

Epoch (phase)" Flux [mJy]* a'^^ 

4.8 GHz 8.6 GHz 



l(-7) 


1.02 (0.17) 


0.82(0.18) 


-0-38!?:?^ 


9(-4)'' 


7.07 (0.16) 


0.57 fO.75) 


7 77-0.49 
-1.93 


2(-3) 


0.75 (0.14) 


< 0.48 (0.16) 


-0.88+"" 


10(+0) 


0.70 (0.14) 


0.72 (0.13) 


+0 05+°-^2 


3(-Hl)'' 


(0.18) 


(0.15) 




7(+ir 


0.67(0.13) 


0.97 (0.14) 


+0 5?+' '"* 


8(H-16/ 


0.51 (0.09) 


0.63 (0.09) 


+0 36^°-'^^ 



° Phase indicates "days from periastron". 

' Numbers in parenthesis indicate the rms value in the maps. 
Spectral index, a = log(F8 f,/F4.8)/ log(v8.6/v4.8) 
Orbital phase also observed in epoch 6 with a compact 
configuration, and hence not considered useful for this project. 

^ Orbital phase also observed in epoch 3 without detection. See text. 
Note that for this epoch Stokes Q at 8.6 GHz is 0.38 mJy, while 
Stokes V and U were not detected. No polarisation signal was 
detected at 4.8 GHz. 

^ Maximum and minimum values of the spectral index when 
adding/subtracting Icr rms to the listed flux. 
As mentioned in the text, we do not detect BP Cru at this epoch. 



Table 4: Characteristics of the objects detected around BP Cru. 
Coordinates and fluxes refer to epoch 8. 





RA(J2000) 


DEC(J2000) 


Flux 


Remarks 








4.8/8.6 [mJy] 




A 


12:26:23.9 


-62:45:58.75 


1.53/1.10 


2MASS'' 


B 


12:26:32.9 


-62:44:55:59 


1.16/0.96 


2MASS'' 


C 


12:26:28.6 


-62:46:17.58 


0.70/0.67 


GSC2.3'' 



" within 7.5" , * within 9.5" , within 1.2" 



(1934-638 for epochs 3 and 4). The flux scale was set by the 
Stokes I fluxes at 4.8 and 8.6 GHz of the phase calibrator 1 148- 
671, i.e. 1.43 and 0.92 Jy, respectively. All solutions obtained to- 
ward the calibrators were applied to the target and images were 
made using different restoring beam sizes to reflect the differ- 
ent array configurations of ATCA in the different epochs. Our 
short observations gave elongated beam sizes, the average beam 
sizes for the epochs with detections (see Table[3]) were 15 x 1.5 
and 10 X 0.8arcsec at 4.8 and 8.6 GHz, respectively. The inte- 
gration time on target varied between 36 and 153 minutes (see 
Table |2]i, which allowed us to reach rms levels in the final maps 
between 0.16 and 0.09 mJy beam"' at 4.8 GHz and between 0.17 
and 0.09mJy beam ' at 8.6 GHz. Even epoch 8, our longest ob- 
servation of BP Cru, was too short and covered a too short range 
of hour angles to allow a safe correction for polarisation leakage. 
What we report here in terms of polarisation is a crude estimate 
done for only one epoch. 



3. Results 

A summary of the detections of BP Cru is presented in Table [3] 
BP Cru was clearly detected at 4 epochs, corresponding to 4 dif- 
ferent orbital phases: two before periastron, one after periastron 
and one at periastron. The clearest detection was obtained in our 
longest integration on-source at epoch 8, 16 days after perias- 
tron. The cleaned map at 8.6 GHz is presented in Fig. 2. All other 
detections are significant to at least 5crrms. We report also on 
two more possible detections, at epochs 7 and 9. The uncertainty 
of the detections comes from the fact that the same orbital phase 
was observed twice, with a detection in only one of the two oc- 
casions. Epochs 6 and 9 correspond to 4 days before periastron, 
but only epoch 9 was useful (see Sect.|2]i; epochs 3 and 7 probed 
BP Cru one day after periastron but a detection could be seen 
only in epoch 7. Apart from these uncertainties, the most strik- 
ing aspect here is that BP Cru seems to emit in the radio regime 
regardless of the orbital phase and especially when it was least 
expected, i.e. close to apastron. 

Another important observation from Tab.[3]is that a appears 
to change from being clearly negative right before periastron to 
being weakly positive after periastron. As visible from the mini- 
mum and maximum values listed in Table[3]Q' seems to be plau- 
sibly positive at and after periastron, while it is possibly negative 
before that. Notice that BP Cru was not detected at 8.6 GHz in 
epoch 2 and we therefore consider the 3crrms level as an upper 
limit for the flux. Furthermore it is important to notice that the 
epochs with possibly negative spectral index coiTelate with the 
X-ray outburst, w hich in its turn is k nown to occur some days 
before periastron (IWatson et al.lll982h . Possible explanations of 
this behaviour are discussed in the next section. 

Finally, from polarisation analysis of the longest observation 
of BP Cru we notice that its radio emission is probably linearly 
polarised, as Stokes Q and Stokes V=U=0. This result is 
to be taken with caution because our observation at that time 
was shorter than three hours, barely enough to correct for polar- 
isation leakage between the ATCA orthogonal feeds. Dedicated 
observations are necessary to retrieve more reliable information 
about the polarisation of the radio emission from BP Cru. 

Several other sources were detected in the ATCA field, these 
are listed in Table |4] Sources A and B where detected at all 
useful epochs. Source A appears to lie some 7.5 arcsec away 
from a 2MASS catalogued object, 12262412-6246060, with 
7//7ir magnitudes 15.83, 15.02, 18.3, respectively. Source B lies 
some 9.5 arcsec from a 2MASS catalogued source, 12263253- 
6244464, with JHK magnitudes 15.67, 14.20, 13.60, respec- 
tively. Source C lies at about 1 .2 arcsec from an object listed 
in the Guide Star Catalogue (GSC)2.3 and has a V magnitude 
of 16.71. The closest 2MASS object to C is 12262821-6246162 
at 2.8 arcsec distance, and having ///TiT magnitudes 15.03, 14.64, 
14.42, respectively. Taking into account the orientation of the ar- 
ray, and considering the large beam sizes in our observations, our 
positions of sources A, B and C are consistent with those of the 
catalogued objects mentioned above. The flux densities of these 
sources remained constant within a; 15%. A further discussion of 
these sources though goes beyond the scope of this paper. 

4. Discussion 

A spherical stellar wind from Wray 977 is expected to gener- 
ate continuum fluxes of 1.9 and 2.6 mJy at 4.8 and 8.6 GHz, 
res pectively, assuming a distance of 3 kpc f or Wray 977 
(P anada & Fellil Il975fc IWright & Bariowl 119751: IScuderi et all 
11991 7 which is more than a factor of two higher than what we 



4 



M. Pestalozzi et al.: Radio emission from the high-mass X-ray binary BP Cru 



-62 44 45 



45 00 



o 

5 



3 
u 




12 26 40 38 36 34 32 30 28 26 24 
RIGHT ASCENSION (J2000) 



flux of its transient radio component is at least 200 times fainter 
than that of GRS 1915+105. 

A problem with this interpretation is that the detection at 
epoch 1 occurs so close to the onset of the X-ray outburst that 
one might question whether a jet can have formed already at that 
time. An alternative source of the non-thermal electrons would 
be a shock that forms in the stellar wind as the neutron star ap- 
proaches Wray 977 and distorts the stellar wind through its grav- 
itational field. 

In the future one should follow the evolution of the radio 
flux from BP Cru over all orbital phases, and every observation 
should be sufficiently long to allow the study of the polarisation 
of the radio emission. This is readily possible using the Compact 
Array Broadband Backend (CABB), newly installed at ATCA. 
Its wide band increases the overall sensitivity by a factor 4, en- 
suring the detection of BP Cru at any orbital phase. Also, with 
sufficiently long integrations to give good signal-to-noise ratios, 
it will be possible to obtain a full spectrum in two bands with 
ideaUy one point every 250 MHz. 



Fig. 2: Stokes I maps of the data taken in epoch 8, at 8.6 GHz. 
The beam size for this epoch was 3.4 x 0.7 arcsec. The box indi- 
cates BP Cru, the circles indicate sources A (to the right), B (top 
in the figure) and C (lowest in the figure) as in Table |4] 



detect. There are several ways in which the predicted fluxes can 
be reconciled with the observed fluxes. Firstly, there are large 
uncertainties in the distance to BP Cru, and a distance of 4- 
5 kpc, sufficient to decrease our observed fluxes with a factor 
of two , is still perfectly compatible with the most recent esti- 
mates (iKaper et al.ll2006l) . Secondly it is not clear that Wray 977 
is emitting a sufficient number of Ly man photons to keep the 
wind completely ionised. In fact, using IScuderi et all (1 19981) we 
estimate that Wray 977 produces some 3-5 times less photons 
than what is required to keep the entire wind ionised. Also, it is 
to be noticed that the radio flux from a stellar wind can vary by 
a factor two (IScuderi et al.lll998h . 

Typically, the radio emission from t he stellar wind of an early 
type star has a positive spectral index jScuderi et alJI 19981) . We 
suggest therefore that the radio emission at epochs 7 and 8 is 
due to the stellar wind from Wray 977. Consequently, the in- 
creased radio flux at 4.8 GHz at the time of the X-ray outburst is 
probably due to a transient, non-thermal component with a nega- 
tive spectral index. There is not enough data, to tell whether this 
component is the most strongly correlated with the orbital phase 
or the timing of the X-ray outburst. 

We propose that the transient non-thermal component ap- 
pearing during the X-ray outburst originates from a weak and 
short lived (relativistic?) jet, as in a microquasar. To test this 
hypothesis, we isolated the transient component by subtract- 
ing the flux from epoch 8 from that at epoch 1. The transient 
component has then fluxes of 0.51 and 0.19mJy at 4.8 and 
8.6 GHz, and a spectral index a = -1.7;^^q. This suggests a 
somewhat steeper index than is usually seen in microquasars 
(e.g. McClintock et al. 2009; Trushkin et al. 2008, in which a > 
-1, but the uncertainties are so large that no firm conclusions 
can be drawn. We also compare BP Cru with the microquasar 
GRS 1915-H105. During an extended period in 2008 Swift/BAT 
registered an X-ray flux ofO.lph cm"- s"' from GRS 1915+105, 
and at th e same time its radio flux at 4.8 GHz remained above 
100 mJy dTrushkin et al.ll2008l) . The X-ray flux from BP Cru at 
the peak of an outburst is even higher, 0.2 phcm"^ s ', but the 



5. Summary and Conclusions 

We summarise our findings as follows: 

- We have presented the first detection of the high-mass X- 
ray binary BP Cru in centimetre radio continuum at 4.8 
and 8.6 GHz. The source has been detected several times at 
different orbital phases, indicating that BP Cru emits radio 
waves regardless of the orbital phase; 

- The radio emission from BP Cru shows a varying spectral 
index, suggesting that it is the superposition of two compo- 
nents: a persistent thermal emission from the wind of the B 
hypergiant mass donor Wray 977, and an episodic emission, 
perhaps a weak jet, that appears at the time of the X-ray out- 
burst; 

- Further observations using the newly installed Compact 
Array Broadband Backend (CABB) will greatly improve 
both the estimate of the spectral index and the information 
about the polarisation of the radio emission of BP Cru. 
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